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1 Introduction 



The direct detection of galactic dark matter (DM) particles in the laboratory using shielded 
nuclear recoil detectors is among the most challenging and high priority goals of astroparticle 
physics. Presently several experiments are taking data with the primary aim of securing a 
convincing detection and thereby confirming the hypothesis that DM is in fact made of relic 
particles. Several ambitious next-generation experiments are planned and from these we hope 
to be able to extract the properties of the DM particle, most notably its mass and cross- 
section for interactions with nucleons. However, there are a number of uncertainties in the 
interpretation of such experimental data. Foremost among these is our poor knowledge of the 
distribution and dynamics of DM in our Galaxy. In particular, the DM velocity distribution 
f{v) strongly affects the conclusions we can draw from the data, for example whether the 
results of different experiments are mutually consistent [1-9]. 

Results from direct detection experiments are usually presented as a signal or exclusion 
curve in the parameter plane of DM mass versus its scattering cross-section, assuming a 
Maxwell-Boltzmann (M-B) velocity distribution with a cut-off at the escape velocity from the 
Galaxy. This is also known as the Standard Halo Model (SHM). There are more sophisticated 
dynamical models of halo DM but to study the effects of these the plot has to be reproduced 
separately for all such velocity distributions which can be cumbersome. 

One might wonder if the experimental data can itself be used to infer the DM velocity 
distribution. Unfortunately, this turns out to be rather difficult [10, 11] because a nuclear 
recoil of energy En can originate from any DM particle which has a velocity larger than some 
minimum velocity VminiE^). Consequently, rather than probing the velocity distribution, di- 
rect detection experiments actually measure the velocity integral: g{vmm) = . f{v)/vd^v. 

Hence, it was suggested [12] that results from one experiment be converted into fmin- 
space in order to predict the event rate in a second experiment, without having to make any 
assumptions concerning the astrophysics. While this approach works well in comparing two 
experiments that are sensitive to a similar range of DM velocities, it is difficult to compare 
experiments that probe different regions of the velocity distribution. It would therefore be 
preferable to have a framework in which the astrophysical unknowns are made explicit so one 
can assess their impact on the extraction of DM properties. 

In this paper we investigate in detail the idea, also suggested in Ref.[12], of mapping 
all experimental results into VmmSpace. To do so, we must know the properties of the DM 
particles — we will assume therefore that their mass can be inferred either from collider 
experiments or the combined information from several direct detection experiments (see e.g. 
Ref.[10, 13]). In this sense, our approach is complementary to the usual analysis of direct 
detection experiments, wherein the astrophysical parameters are held fixed and the DM mass 
and cross-section are allowed to vary. 

Just as the usual presentation of results shows whether all experiments are consistent for 
some range of DM mass and scattering cross-section, our treatment will illustrate whether it 
is possible to find a DM velocity distribution that can reconcile all such experimental results. 
To demonstrate this, we focus on the DM mass range 6-15 GeV motivated by recent claims 
of signals from DAMA, CoGeNT and CRESST-II [14-16] and show that (7(fmin) can indeed 
be chosen so as to enable a consistent description of these experiments. However, upper 
limits on signals from XENON [17] and CDMS [18, 19] cannot be evaded by appealing to 
astrophysical uncertainties alone. If all these experimental results are correct then to resolve 
the tension between the signals and upper limits requires non-standard DM interactions. 
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This paper is organised as follows: Section 2 introduces notation and sets out the 
analysis framework. In Section 3 we discuss how a measurement of the differential event rate 
in a direct detection experiment can be translated into a measurement of g{vmin) and how 
an experiment that does not observe a DM signal can constrain the velocity integral. We 
present measurements and limits on g(wmin) for the most relevant direct detection experiments 
in Section 4 and also assess the compatibility of CoGeNT and CRESST-II. In Section 5 we 
do a similar analysis for experiments which observe an annual modulation of the event rate 
and assess the compatibility of DAMA and CRESST-II. In Section 6 we discuss briefly 
how varying the particle physics properties of DM can affect the analysis, and present our 
conclusions in Section 7. 



2 The velocity integral 



The differential event rate in the laboratory frame for the scattering of DM particles on nuclei 
is given by 

-CUA,Z)F\EMvmin) , (2.1) 



d^R 2m^^2^ 

where p is the DM density, is the DM mass and /Xnx is the reduced DM-nucleon mass. 
The 'velocity integral' g(wmin) is defined by 

,K.„.*).r^i^±^d3., (2.2) 

where f{v) is the local DM velocity distribution evaluated in the galactic rest frame, v = \v\ 
and v-p^{t) is the velocity of the Earth relative to the galactic rest frame [20, 21] . The minimum 
velocity required for a DM particle to transfer an energy E-^ to a nucleus is 




{Ek) = x^-^, (2.3) 



where m-N is the mass of the target nucleus and jj, the corresponding reduced mass of the 
DM-nucleus system. 

To avoid clutter in our formulae, we define Ct(^, Z) = {fp/ fnZ + {A — Z)), where A 
and Z are the mass and charge numbers of the target nucleus and fn,p denote the effective 
DM coupling to neutrons and protons, respectively. We assume fn/ fp = 1 unless explicitly 
stated otherwise. Finally, o"„ is the DM-neutron cross-section at zero momentum transfer 
and F{En) is the nuclear form factor, which encodes the loss of coherence as the momentum 
transfer deviates from zero. We use the Helm form factor from Ref.[22]. 

For our purposes, it will be convenient to absorb the DM mass, cross-section and density 
into the definition of the velocity integral and consider the rescaled velocity integral 

g{Vrain) = g{Vmm) , (2.4) 

which has a typical value in the range 10^^^ — 10~^^ day~^ (see Figure 1). 

At this point, it is worth developing some intuition for g{vmin) and elucidating its 
relation to the more familiar DM velocity distribution f{v). In the left panel of Figure 1, 
we show the canonical M-B velocity distribution of the SHM, along with a sharply peaked 
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Figure 1. The one-dimensional velocity distribution 47rw^/(w) (left) and the rescaled velocity integral 
5(i'min) (right) for the Standard Halo Model (blue, dotted) and a sharply peaked DM stream (green, 
solid). To calculate ^(t'min), we have assumed cr„ = 10^''*' cm^, p = 0.3 GcV/cm"^ and — 15 GeV. 
The velocity of the DM stream is Vg — (250, 0, 0) km/s, giving Vmin = \vs — ve\ ^ 330 km/s. 



function which could arise from a stream of DM or alternatively can be thought of as a shifted 
M-B distribution in the limit of small velocity dispersion. In the right panel we show the 
corresponding g{vmm)-^ As observed earlier [23], g{vmm) has the important property that for 
any velocity distribution it is a decreasing function of Vmm since f{v) > 0. As the velocity 
dispersion decreases, g{vmin) becomes more like a step- function which is the minimal form 
consistent with this property. 

3 Measurements and constraints 

We now discuss how information on ^(fmin) can be extracted from direct detection experi- 
ments.^ First we consider experiments that observe a signal that can be interpreted in terms 
of DM scattering, namely DAMA, CoGeNT and CRESST-II. In this case, the measured sig- 
nal can be directly translated into a measurement of ^(fmin)- Then there are experiments 
that do not observe a DM signal above their expected background and thus constrain the 
differential DM scattering rate which converts into a bound on ^(fmin)- Since ^(fmin) may be 
time-dependent in principle, the resulting limit is strictly speaking on the average of ^(vmin) 
over the period the corresponding experiment was taking data. 

3.1 Measuring the velocity integral 

It is relatively straightforward to infer the value of g{vmm) at Vmm{E^) from a nuclear recoil 
detector such as CoGeNT that consists of a single element and measures the differential event 
rate dR/dEn at the recoil energy Sr. Inverting Eq. (2.1), we find 

2 1 d^ 



For a real detector we must take the detector resolution IS.E and the detector efficiency e{E-^) 
into account. If the detector measures the differential event rate in the interval [Ei,E2], we 



^We neglect the time dependence of gi(i;min) here, which would lead to a slight smearing of the step function, 
similar discussion is presented in the appendix of Ref. [12]. 
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can infer the value of ^(fmin) for Umin in the interval [fmin(-E'i — AE),Vram{E2 + ^E)]: 

2 1 dR 



CUA,Z)F^EMEn) dEn ■ 



The recoil energy where the nuclear form factor and the efficiency are evaluated should in 
the range Ei < E^ < E2; we assume that the energy bin is sufficiently small that these 
quantities do not vary significantly with Eji within the bin.^ 

Additional care is required for experiments like DAMA and CRESST-II which have 
targets with more than one element. The DAMA detector has both sodium and iodine 
nuclei, the latter being much heavier than the former (mj ~ 5.5mNa)- For the light DM we 
consider, the kinetic energy of the recoiling iodine nucleus is below the energy threshold of 
the detector so it is a good approximation to assume that all contributions to dR/dE^^ come 
from scattering on Na nuclei and use Eq. (3.2) taking the mass fraction of Na in the detector 
into account (see Section 5). 

For CRESST-II which contains oxygen, calcium and tungsten, the situation is more 
complicated. Since tungsten is much heavier than oxygen and calcium, we can ignore it as 
the corresponding recoil energy is below the detector's energy threshold. However oxygen 
and calcium are relatively close in mass so extracting information on ^(t'min) becomes more 
complicated. We discuss these difficulties and how to deal with them in Appendix A. 

3.2 Constraining the velocity integral 

Next we consider experiments that do not measure the differential event rate but are rather 
able to bound dR/dE^ over some energy range. This is usually converted into a bound on 
an for a given DM mass, assuming a galactic halo model such as the SHM. Instead we wish 
to use the bound on dR/dE^i to constrain ^(fmin) for different values of Umin-^ 
Since f{v) > 0, the following inequality holds for any value of Vmin- 

g{Vmin) > g{v^mn) ©("t^min " I'min) , (3.3) 

where ^({'min) is a constant and Q{x) is the Heaviside step function. Of all velocity integrals 
that have ^(^min) = 50i the one defined by g{vm\n) = 9o0(^min — ^min) thus predicts the 
smallest number of events in any given experiment. Of course (7(fmin) ©(^min — I'min) is 
not a realistic model for the galactic halo, nevertheless it is a valid velocity integral that can 
be used to predict the event rate for a given experiment. A more realistic halo must satisfy 
Eq. (3.3) and will therefore necessarily predict a larger event rate. Consequently, if we can 
reject the case when ^(fmin) = 9o0('^'min — Vmm)-, we can also reject any other other halo 
model giving g{v,ain) = go- 

Hence if an experiment places an upper bound on di?/di?R, we can correspondingly 
bound g{vram) by fixing Vmm = Vmin and finding the smallest value of ^(t'min) such that the 
predicted event rate for g{vram) = 5(^'min) ©('Wmin — ^^min) is larger than the measured value (at 
a given confidence level). Repeating this procedure for all Wmim we obtain a continuous upper 
bound on ^(umin) — see Figure 2. If the exclusion curve thus obtained touches (or crosses) 

^This assumption is good for CoGeNT and DAMA for the light DM that we consider. However it is not 
a good approximation for CRESST-II where e{E^) can change rapidly within a single bin that contains the 
onset of a new detector module; to avoid this we bin the data so such regions are avoided (see Appendix D). 

''It is straightforward to put a bound on g(«min) from experiments that use a single element; difficulties 
arising in experiments with more than one type of nuclei are discussed in Appendix A. 
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Figure 2. An illustration of the method developed in Section 3.2. The blue dashed lines correspond 
to velocity integrals of the form go 6({'min ~ Vmin) for different choices of and 'Omin- The predicted 
event rate in the XENONIOO detector for each of these velocity integrals can be rejected by the data 
from XENONIOO. The green line shows the resulting bound on g(wmin)- 



the g{vmin) curve for a particular f{v), this model will be excluded by the experimental data 
at the same confidence level used to construct the exclusion bound. In other words, if the 
bound on g{vmin) lies below (some of) the values of ^(umin) implied by the measured recoil 
spectrum at DAMA, CoGeNT or CRESST-II, it will not be possible to find a halo model 
that consistently describes all measurements and evades all experimental bounds. 

Note that the converse of this statement is not necessarily true: even if g(vmin) stays 
below a given exclusion curve for all values of Vmm, this does not imply that the model is 
not excluded by the corresponding experiment. By construction, we have ensured only that 
the predicted differential event rate in each individual energy bin is below the respective 
experimental limit. However, if the measured velocity integral remains close to the exclusion 
limit over a wide range of Vmin, the predicted total event rate, i.e. the integral of the differential 
event rate over several bins, may be excluded at the chosen confidence level. 

4 Results 

Now we use the framework developed in the previous section to analyse recent results from 
direct detection experiments (discussed in Appendix D). We use the events reported by 
CoGeNT and CRESST-II to measure giv^ain) and the null results from XENON, CDMS, 
SIMPLE and the CRESST-II commissioning run to constrain it. Our results are shown in 
Figure 3 for various values of the DM mass . As we increase the exclusion curves as well 
as the CoGeNT data points move according to fmin ~ (which holds when < ttin) — 

the CRESST-II data points and exclusion curves scale differently because the mass of oxygen 
is comparable to in the lower panels. It should be clear from this consideration that the 
CoGeNT and CRESST-II data points are always excluded for 6 GeV < < 15 GeV. 

For all masses considered, one can see the obvious conflict between the measurements 
by CoGeNT and CRESST-II and the exclusion bounds from other experiments. In each case, 
most of the measured points lie above one (or more) of the exclusion curves. We conclude 
that a consistent description of all experiments is not possible for any model of the DM 
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Figure 3. Measurements and exclusion bounds of the velocity integral ^(wmin) for different DM 
masses m^. The DM interaction is spin- independent and elastic with /„//p — 1. Values of g(wmin) 
above the lines are excluded with at least 90% confidence. For = 6, 9 and 12 GeV the data points 
for CRESST-II have been obtained using two different methods, as described in Appendix A. It is not 
possible to find any model for the DM halo that provides a consistent description of all experiments. 



halo if the DM particles undergo elastic spin-independent scattering with fn/ fp = 1- There 
is no functional form for g{vrmn) that would allow a DM interpretation of the CoGeNT or 
CRESST-II data consistent with other experiments. 

4.1 A consistent description of CoGeNT and CRESST-II 

Even though we cannot find a halo model that provides a consistent description of all ex- 
periments, it is seen from Figure 3 that CRESST-II and CoGeNT probe ^(I'min) at different 
ranges of fmin- Therefore it should be possible to choose ^(fmin) such that we obtain a con- 
sistent description for these two experiments. This choice must be different from the SHM 
for which the best-fit DM regions of CRESST-II and CoGeNT do not overlap [16, 24]. 

Of course, we cannot vary ^(t'min) arbitrarily — in the end, the velocity integral must 
arise from a reasonable self-consistent model of the DM halo. Therefore in Figure 4, we 
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Figure 4. The velocity integral ^(wmin) for the SHM with different parameters (left) and alternative 
descriptions of the DM halo introduced in Appendix C (right). The shaded region corresponds to the 
values of (7(wniin) observed in the GHALOs simulation [5]. 



examine the range of predictions for ^(f min) from a variety of reasonable models of the galactic 
halo. First we note that g{v^^^) can change considerably even in the context of the SHM, if we 
vary vq and v^sc within their observational bounds (see the left panel of Figure 4). Secondly, 
the SHM is unlikely to be an accurate description; indeed many alternative models and 
parameterisations for the halo exist in the literature. We present an overview in Appendix C 
and show the corresponding velocity integrals in the right panel of Figure 4. Even for a 
fixed choice of v^sc and vq there is a large spread in the predictions of the velocity integral, 
especially close to the cut-off. 

We conclude that given the spread in the predictions for ^(umin) it should be possible 
to bring CoGeNT and CRESST-H into better agreement. In this context it is instructive to 
ask why CRESST-II favours larger DM masses than CoGeNT (see left panel of Figure 6). 
To understand this, we show the SHM prediction for ^(^min) in Figure 5 together with 
the measurements from CoGeNT and CRESST-II for m 



X 



9 GeV. We observe that the 
SHM prediction of the rescaled velocity integral below 500 km/s is slightly too flat to fit the 
CoGeNT data, thus favouring smaller DM masses. Moreover there is an additional constraint 
(which is not apparent in Figure 5) from the fact that CoGeNT does not observe a signal at 
higher energies. This constraint additionally disfavours the case m-^ > 9 GeV if the rescaled 
velocity integral is too flat. Consequently if we want to push the CoGeNT region to larger 
DM mass, we need to make the rescaled velocity integral steeper below 500 km/s. As a 
result our new velocity integral will predict fewer events from oxygen in the lowest bins of 
CRESST-II because ^(umm) becomes smaller around 500 km/s. To compensate, we must 
increase the contribution from calcium in the lowest bins, corresponding to larger values of 
^(^'min) around 600 km/s (see also the discussion in Appendix A). 

We can easily achieve this goal by adding two SHM-like contributions — one with 
relatively small velocity dispersion and low cut-off which dominates at low values of Vmin, 
and a second one with much larger vq and Vesc- Of course, the resulting velocity distribution 
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Figure 5. Measurements of the velocity integral ^(wmin) compared to the predictions of the SHM and 
the Modified Halo Model. In both cases, we have assumed (t„ = 10"''^° cm^ and p = 0.3GeV/cm'^. 
The Modified Halo Model provides a consistent description of the data from CoGeNT and CRESST-H 
for a DM particle with — 9 GeV. 



is not a self-consistent model of the DM halo; it is only a convenient parameterisation to 
illustrate how the velocity integral can be changed to simultaneously accommodate several 
direct detection experiments. Nevertheless, it is not inconceivable that such a velocity integral 
can be realised in models with a strong anisotropy, where the tangential component with low 
velocity dispersion dominates at low energies and the radial component with high dispersion 
dominates at high energies. We shall refer to this new velocity integral as the Modified Halo 
Model (MHM). 

We show the corresponding best-fit regions and exclusion limits in the - an plane 
in Figure 6. As desired, the CoGeNT and CRESST-II region^ now overlap at ~ 9 GeV. 
The common parameter region is, nevertheless, clearly excluded by XENON and CDMS, an 
observation which we could equally well have made from Figure 3. 

5 Measurements and constraints for modulation amplitudes 

We have seen that for all the cases considered, there is strong tension between the exclusion 
limits from null results and the event rate observed by CoGeNT and CRESST-II. This seems 
to suggest that at least some of the events seen by CoGeNT and CRESST-II may not be 
due to DM scattering but a new source of background [25]. We will therefore now consider a 
more specific signature of DM interactions, namely the annual modulation of the signal due 
to the motion of the Earth relative to the galactic rest frame. This will help to identify the 
DM signal if the background is not time-dependent. 

In fact, two experiments have observed an annual modulation of their signal, namely 
DAMA [14] and CoGeNT [15]. By considering the modulation in fmin space, it was observed 
[26] that the two modulation signals are (marginally) compatible with each other, independent 
of astrophysics. Now we wish to include these annual modulations into our analysis of 

^For CRESST-II we show the la best-fit region using the publicly available data as stated in Appendix D. 
Note that this corresponds roughly to the 2a region published by the CRESST collaboration [16]. 
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Figure 6. 90% confidence regions inferred from tlie absolute rate observed by CoGeNT and la 
regions from CRESST-II as well as the 90% exclusion limits from other experiments for the SHM 
(left) and the Modified Halo Model (right). By construction, the Modified Halo Model leads to an 
overlap of the CoGeNT and CRESST-II regions. 



5(^min) to determine whether these signals are compatible with the null results from other 
experiments. We first discuss these constraints without making any assumptions about the 
modulation fraction and then discuss how the modulation fraction can itself be reasonably 
constrained in order to obtain stronger experimental bounds. 

The time dependence of g{vra\n-,t) can be approximately parameterised as 



g{Vmm,t) = g{Vmm) 



1 -1- A{Vmm) ■ COS 271- 



1 yr 



(5.1) 



where (7(fmin) is the time average of g{vram,t)- Note that as emphasised earlier [5], both the 
modulation fraction A and the peak date to can depend in general on Wmin- For future use 
we define the modulation amplitude Ag(umin) as 

Ag{Vmm) = ^ [g{Vmm, to{Vmm)) " ^(^^mm, io(^^min) + 0.5 yr)] = A{Vjnm)g{Vmin) , (5.2) 

and introduce the rescaled modulation amplitude: Ag{vmin) = 'Ag{vmm)(^nP/fnx- 

First we need to constrain the function to(^min)- CoGeNT and DAMA probe largely 
the same region of Vmm space, so it would be inconsistent to assume different phases for the 
two experiments. It was shown [26] that the phase measured by DAMA and by CoGeNT 
are consistent at the 90% confidence level for sufficiently low DM mass. We will there- 
fore take to{vram) = 146 days which is the best-fit value from DAMA.^ As the modulation 
fraction satisfies A < 1, we see from Eq. (5.2) that A^(t;min) < g{vmm)- We can easily check 
whether DAMA and CoGeNT satisfy this inequality by plotting the respective measurements 
of A^(fniin) on top of the constraints for ^(i^min) from other experiments which we obtained 



^Varying to within the la region allowed by DAMA changes the best-fit value for the CoGeNT modulation 
amplitude by about 20%, which does not significantly afi^ect any of our conclusions. 





Figure 7. Measured values of Ag{vyain) from DAM A and CoGeNT compared to the exclusion limits 
from other experiments. For the upper panels, no assumptions on the modulation fr'action have 
been made, for the lower panels, we assume that the modulation fraction is bounded by the red 
line in the right panel of Figure 8. Even for weak assumptions on the modulation fraction, there 
is significant tension between the different experiments, most notably it is impossible to find a DM 
velocity distribution that describes the observed modulations and evades the bound from XENONIOO. 



above. In complete analogy to Section 3.1 we can calculate Ag{v^i^) from the modulation 
of the event rate, AdR/dE^i, seen by DAMA and CoGeNT 

1 d-R 

A^(._) = 2^1^ cUA,Z)F^iE^) ^dE^ ' ^''^^ 

The resulting values of Ag^Vmin) are shown in Figure 7. We observe that the bounds from 
the XENONIO and XENONIOO experiments remain strong even if we do not make any as- 
sumptions about the modulation fraction. However other experiments do not significantly 
constrain Ag{v^in)- We consider therefore whether it is reasonable to make stronger assump- 
tions about the modulation fraction and thus obtain more stringent experimental bounds. 
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Figure 8. Left: The velocity integral and the modulation fraction for different descriptions of the 
galactic DM halo. Right: The constraint for the modulation fraction extracted from the CoGeNT 
data (boxes) together with the bound that we will adopt (red solid line). The shaded region in both 
panels corresponds to the modulation fractions observed in the GHALOg simulation [5] . 



5.1 Constraining the modulation fraction 

We will now discuss what can be reasonably assumed about the modulation fraction given 
known models of the galactic halo, and how it can be constrained once the velocity integral 
has been measured. The predicted modulation fraction for various halo models are shown in 
the left panel of Figure 8. We observe that for most values of fmin it is significantly below 
100%. Note that a modulation fraction of 100% implies that no signal is observed at to + 0.5 
yr, which is possible only if Vmin > Vesc + ^eC^o + 0.5 yr). 

The empirical bound that we will adopt is shown as the red line in the right panel of 
Figure 8. Note that this is a weak bound in the sense that we allow modulation fractions 
much larger than the value expected in any reasonable description of the DM halo. If we 
require that the modulation fraction satisfies A{vmm) < ^"^'^^(^^min)) any direct detection 
experiment that constrains g{vrain) < ff™'^^(^'min) will also constrain Ag{vram) according to: 

^giVmin) = A{Vn,in)g{Vmm) < ^'^''''(^^mm)5"'^''(^'min) • (5.4) 

The resulting measurements and constraints for A^(?;min) are shown in the lower panels of 
Figure 7. Even with our weak assumption concerning the modulation fraction, the bounds 
on A^(ujnin) have become significantly stronger. Decreasing the mass of the DM particle now 
helps significantly because at high values of Vmin larger modulation fractions are allowed so 
the exclusion limits become less stringent. However for small masses it becomes increasingly 
difficult to explain the modulation observed by CoGeNT above 2 keVee because these recoil 
energies correspond to very large velocities, where we expect ^(fmin) and therefore Ag{vmin) 
to be small. 

We now consider whether the modulation fraction can be constrained further using 
measurements of g{vmin)- Comparing the left panel of Figure 8 to the right panel of Fig- 
ure 4, we observe that there is a direct correspondence between the velocity integral and the 
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Figure 9. Measurements of the modulation amplitude Ag{v,nin) compared to the predictions of 
the SHM and the Modified Halo Model. In both cases, we have assumed (t„ — 10^'*" cm^ and 
p — 0.3 GeV/cm^. While in the SHM, the predicted modulation amplitude for these parameters is 
too small to explain the observed modulations, we can get a sufficiently large modulation amplitude 
in the Modified Halo Model. 



modulation fraction: the steeper the slope of the velocity integral, the larger the modulation 
fraction. This result is intuitively clear — if ^(fmin) changes rapidly with varying v^iin we 
can also expect large changes between the rates observed in summer and winter. This can 
in fact be used to constrain the modulation fraction once the velocity integral is known. We 
derive in Appendix B that 

^ _ gjvmm, summer) - g{vmm, winter) ^ gjvmm - u) - gjVrmn + u) 
5(Vmin, summer) +5 (i-min, wmter) 2g{Vmm) 

where u = 29.8 km/s if we assume an anisotropic halo and u = 15.0 km/s if we assume an 
isotropic halo. 

As an example, we can use the CoGeNT data to derive a bound on the modulation 
fraction for a given DM mass. For this purpose, we need to bin the measured events in 
such a way that the resulting bin width in Wmin-space is equal to u (see Appendix D). The 
resulting constraints are shown in the right panel of Figure 8. We observe that assuming an 
isotropic halo, these constraints are actually quite severe, limiting the modulation fraction 
to 20% at 

Vmin — 400 km/s, which in fact coincides with our adopted bound for A{v^\Yi). For 
an anisotropic halo, modulation fractions of up to 40% may be possible in the same range. 

This observation is also relevant for a self-consistent DM interpretation of the CoGeNT 
experiment. To bring the modulation observed by CoGeNT into agreement with the measured 
event rate, a relatively large modulation fraction is required [27]. If CoGeNT were to reject a 
significant part of its signal as background, as suggested in Ref.[25], even larger modulation 
fractions will be necessary for self-consistency. The constraint on the modulation fraction 
derived above can then be used to conclude that the DM halo must be anisotropic. 

5.2 A consistent description of DAMA and CRESST-II 

As can be seen from Figure 7, DAMA and CoGeNT probe the same region of ?;min space, 
so it is not possible to improve their agreement by changing the halo model as we did for 
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[GeV] [GeV] 

Figure 10. 90% confidence regions from CoGcNT (inferred from both tlie modulation and the 
absolute rate) as well as the la region from CRESST-II and the DAMA region for the SHM (left) and 
the Modified Halo Model (right). The 90% exclusion limits from other experiments are also shown. 
By construction of the Modified Halo Model, the four best-fit regions are now in good agreement. 



CoGeNT and CRESST-II. Fortunately, such a modification is not necessary because DAMA 
and CoGeNT favour roughly the same modulation amplitude. The more interesting question 
is whether this modulation amplitude is compatible with the average value of ^(I'min) inferred 
from CRESST-II and the CoGeNT unmodulated event rate. 

For a cross-section of (T„ w 10"^*^ cm^, the SHM prediction of the modulation amplitude 
is too small to account for the modulation seen by DAMA and CoGeNT [26-30]. This is 
obvious from plotting the SHM prediction for Ag{vmm) together with the measured modu- 
lation amplitudes (see Figure 9). However, we have already noted in Section 4.1 that the 
SHM must be modified if we want to simultaneously explain CoGeNT and CRESST-II. Now 
we consider if we can bring these experiments into agreement with DAMA (and therefore 
also with the CoGeNT modulation) by additionally allowing a larger modulation fraction. 
Of course, for a given velocity integral the modulation fraction must satisfy Eq. (5.5), so we 
cannot choose arbitrarily large values. 

For this purpose, we assume that (7(fmin) is given by the Modified Halo Model from 
Section 4.1 and that the modulation fraction saturates the bound A"^'^^{vmin) shown in Fig- 
ure 8. By the reasoning in Appendix B these two assumptions are compatible only if the DM 
halo is highly anisotropic. As is demonstrated in Figure 9, we then obtain a sufficiently large 
modulation amplitude to describe the DAMA and CoGeNT modulations. Consequently, the 
four best-fit regions in the traditional (T„ — parameter plane are now in good agreement 
(Figure 10). However these are of course excluded by XENON and CDMS. Because of this 
obvious contradiction we wish to emphasise again that we do not consider present data suf- 
ficient to actually determine the velocity integral or the modulation amplitude. We use it 
only to illustrate how our method can be used to bring future, more reliable, data sets into 
agreement. 



- 14 - 




Figure 11. Measurements and exclusion bounds of tlie velocity integral ^(wmin) (left) and the mod- 
ulation amplitude Ag{v^in) (right) under the assumption fn/fp = —0.7. As above, we have assumed 
that the modulation fraction is bounded by the red line in the right panel of Figures 8. The bounds 
from liquid xenon experiments can mostly be evaded in this case but the bounds from CDMS-Si, 
SIMPLE, and the CRESST-II commissioning run become much more important. 



5.3 Additional contributions to the local dark matter density 

Our discussion of g{vmin) is quite independent of the origin of the local DM density. However 
to predict the velocity integral (Figure 4) and the modulation amplitude (Figure 8) we have 
assumed that the local DM density is completely dominated by the contribution from the 
galactic DM halo. In general there may be other significant contributions to the local DM 
density, e.g. from DM 'streams' and a 'dark disk'. We will now briefly discuss how these can 
alter the theoretical predictions of givrain) and A{vrain)- 

Streams 

A^-body simulations show that the DM velocity distribution is not a smooth function but 
instead has a significant amount of sub-structure [31] due to the presence of tidal streams. 
Such streams can have large velocities relative to the local standard of rest and can therefore 
contribute to g{vm\n) at large values of fmin- The result is an edge in the velocity integral (see 
Figure 1) as well as a significant increase in the modulation fraction [24, 32]. Nevertheless, 
the bound on the modulation fraction derived in Appendix B remains valid. In other words, a 
large modulation fraction is always visible as a steep decrease in the velocity integral and can 
therefore be probed by measuring the differential event rate (see also Ref. [33] for a discussion 
on how to constrain DM streams with CoGeNT). 

Dark disk 

Most A^-body simulations of the galactic DM halo neglect the effect of baryons. Baryonic 
matter is expected to effect the formation of a DM disk that can contribute a similar amount 
to the local DM density as the halo [34, 35]. Such a dark disk is expected to have a much 
smaller velocity dispersion than the galactic halo and to be co-rotating with only a small lag. 



- 15 - 




Figure 12. 90% confidence regions from CoGcNT (inferred from both tlie modulation and the 
absolute rate) as well as the la region from CRESST-II and the DAMA region for the SHM (left) and 
the Modified Halo Model (right) for the assumption /„//p = —0.7. The 90% exclusion limits from 
other experiments are also shown. For the MHM, the region of overlap from DAMA, CRESST-II and 
the CoGeNT modulation is only just excluded by current experimental limits. 



Consequently, a dark disk will contribute to g{vram) only for Vmin < 250 km/s — a region that 
is irrelevant in the context of light DM [4, 6, 36]. For DM particles with rriy. > 50 GeV, one 
can however expect a significant increase in both the velocity integral and the modulation 
fraction at low energies. 

6 Varying particle physics 

All the results shown so far correspond to standard assumptions for the interaction between 
DM particles and detector nuclei. Various modifications of these assumptions have been 
proposed to reduce the tension between various experiments. A particularly interesting 
possibility is isospin-dependent DM couplings, i.e. fn fp- Choosing the ratio of couplings 
fn/ fp ~ —0.7 strongly suppresses the bounds from all experiments using xenon [26-30, 37-40] , 
which are otherwise the most constraining (see Figure 10). This observation has generated 
much interest in DM models leading to such couplings [40-44]. 

We show the resulting measurements of ^(fmin) and Ag{vmin) in Figure 11. The bounds 
from XENONIO and CDMS-Ge can be evaded in this case but the bounds from CDMS-Si, 
SIMPLE, and the CRESST-II commissioning run become more important. Also, for fn/ fp — 
—0.7 the agreement between CoGeNT and CRESST-II as well as the agreement between 
the DAMA and CoGeNT modulations get considerably worse. Interestingly, the agreement 
between CoGeNT and CRESST-II would improve considerably, if the CoGeNT event rate 
were to be reduced through to the subtraction of an additional background. The same 
observation can also be made from Figure 12 (see also Refs.[24, 45]). For the MHM, the region 
of overlap from DAMA, CRESST-II and the CoGeNT modulation is only just excluded by 
current experimental limits. 
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7 Conclusions 



We have discussed how measurements of the differential event rate in nuclear recoil detectors 
can be used to extract information on the velocity integral ^(t^min)- By converting experi- 
mental results into t;min-space, we can compare the results from various recent DM direct 
detection experiments without making any assumptions concerning the properties of the DM 
halo. This strategy has several direct applications. 

The first is to compare experiments that observe a potential DM signal with the null 
results from other direct detection experiments. To do this we need to translate bounds on 
the differential event rate into bounds on g{vmm)- In order to exclude a certain value g{vmm) 
we must demonstrate that any velocity integral that takes this value necessarily predicts an 
unacceptably large number of events in at least one of the experiments. We have utilised 
this approach to demonstrate that astrophysical uncertainties are not sufficient to solve the 
present conflicts between direct detection experiments. Specifically, the values of ^(fmin) 
favoured by CoGeNT and CRESST-II are excluded by XENON and CDMS. 

The second application is to assess the compatibility of several experiments that observe 
a positive signal, independently of astrophysical assumptions. If the experiments probe the 
same region of Vmm space, but measure contradictory values of g{vynm), we can conclude that 
the two experiments cannot be brought into agreement by changing astrophysics. However, if 
the two experiments probe different regions of Wmin space, either because they employ different 
targets or different thresholds, it will often be possible to modify the velocity integral in such 
a way that both experiments are brought into agreement — provided that a monotonically 
decreasing function fitting all data exists. 

We have demonstrated this possibility by considering the recent signals seen by CoGeNT 
and CRESST-II. For this purpose, we have developed a method to treat targets consisting 
of several different nuclei. While the DM masses inferred from these measurements do not 
agree for the SHM, we can modify the halo model to make the respective best-fit regions 
overlap. To assess how reasonable such a modification of the SHM is, we have studied a large 
variety of different dynamical models for the galactic DM halo and calculated the respective 
predictions for g(vram)- The corresponding spread in predictions for the velocity integral is 
rather large and justifies considering significant departures from the SHM value. 

Similarly we can convert measurements of the modulation of the event rate into mea- 
surements of the modulation amplitude of the velocity integral A(^(';;jii;n)- Again such a 
conversion can be used to check whether several experiments observing an annual modula- 
tion are compatible and to confront them with constraints from null results. We find that 
the annual modulations observed by CoGeNT and DAMA are consistent with each other but 
in conflict with the bound from XENONIOO even if we allow a modulation fraction of 100%. 
By constraining the modulation fraction with the absolute rate we have shown that such a 
large modulation fraction is in fact inconsistent with the CoGeNT data. Nevertheless for 
anisotropic halos it remains possible to have a sufficiently large modulation fraction to allow 
for a consistent description of CoGeNT, DAMA and CRESST-II, although the bounds from 
other experiments then become much stronger. 

In conclusion, for standard interactions between DM particles and nuclei, the tension 
between direct detection experiments cannot be completely resolved by varying astrophysical 
parameters. If the signals reported by DAMA, CoGeNT and CRESST-II are interpreted in 
terms of DM while all null results are taken at face value, non-standard interactions like 
isospin-dependent DM have to be invoked. We have discussed this briefly and found that 
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for fn/ fp = —0.7 experimental constraints are significantly weaker and can be evaded by 
modifying the halo to have a larger modulation fraction. In this case, the CoGeNT and 
DAMA modulations as well as the CRESST-II event rate can be brought into agreement with 
all other experiments except the bound from SIMPLE and the CRESST-II commissioning 
run. We look forward to more data to help resolve the present puzzling situation. 
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Appendix 

A Including different target elements 

As demonstrated in Section 3.1, it is relatively easy to convert energy spectra into Umin-space 
for a target consisting only of a single element. However many experiments, for example 
SIMPLE and CRESST-II, employ a combination of different elements as target. In this case, 
a nuclear recoil of a given energy corresponds to different t'min depending on the recoiling 
nucleus. 

We discuss how to disentangle the arising ambiguities. Note that these considerations 
may also be important for targets consisting of different isotopes of a single element, especially 
when considering isospin-dependent interactions. 

A.l Constraining the velocity integral 

First, we discuss how an upper bound can be placed on the velocity integral for a null result 
from an experiment consisting of more than one element. The most obvious example is the 
SIMPLE experiment, consisting of C2CIF5. We apply the method presented in Section 3.1 
separately to each element (or isotope) by requiring that the number of recoil events expected 
for this element alone does not significantly exceed the total number of events observed. 

For each element we then obtain a bound g^^^{vmiii) on the rescaled velocity integral. 
We can combine these bounds into a total bound according to 



This procedure is illustrated in Figure 13. 

In fact, it is sufficient to calculate 5™ '"(t'min) for only one of the elements in the target, as 
we can use this bound to infer the bounds for all other elements. Suppose we have determined 
the bound g{vyam) < 5™^^(^^min) for ^^min = ^"1 for an isotope with charge Zi, mass number 
Ai and form factor Fi. For a different isotope, we can now infer the value of the rescaled 
velocity integral at the minimal velocity V2, which satisfies Ei{vi) = E2{v2), meaning that 




(A.l) 
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Figure 13. The individual bounds on ^(winin) for the different elements and isotopes as well as the 
combined total bound from the SIMPLE experiment. 

the two velocities correspond to the same maximum recoil energy on the respective nuclei. 
The bound from the second isotope is then given by 



;;max/'., x _ ^V^-ry^ I, )^ vy^ vy^ U ) -max/-, n /a r,N 



ciC^{Zi,Ai)Fi{Ei{vi)f 

C2C^{Z2,A2)F2{E2{V2)) 

where ci^2 are the respective concentrations of the isotopes. 



A.2 Measuring the velocity integral 

We turn now to the more difficult task of inferring information on g(wmin) from a potential 
DM signal in an experiment consisting of several different elements. The obvious application 
is the CRESST-II experiment, which we will use here for illustration. We consider only the 
contribution of calcium and oxygen, because for DM masses of 15 GeV and less, nuclear 
recoils on tungsten have recoil energies well below the threshold of the detector. 

The differential event rate in CRESST-II is given by the sum of the differential event 
rates for the individual elements, i.e. 

dRV°' i^^yo) ^ /^y^") ^^^^ 



To cause a recoil of energy Eji, a DM particle must have a minimum velocity of ?;^^(£'r) for 
scattering off oxygen and vl^^^{E^) for calcium. Consequently, if we measure the differential 
event rate at the energy Eji, we simultaneously probe the rescaled velocity integral at both 
of these velocities. Therefore, we can only extract information on 9{v^^^{Eji)) if we know 
the ratio of the differential event rates for scattering off oxygen and calcium 

jdR/dEn)^^^ _ CUZ^^\A(^^).F(o\E^)'.g{vi^iiE^)) 

(dfi/di?R)(C^) C2(Z(Ca),^{Ca)) . i7(Ca) )2 . g^v^^^ (E^)) ' 

We present two different methods to determine this ratio. The first makes use of the fact 
that the rescaled velocity integral is a decreasing function of Wmin to place an upper bound on 
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the contribution from calcium. The second method requires a special binning of the data in 
such a way that the contribution of calcium can actually be determined from the combined 
information of several bins. Once we have an estimate of the differential event rate from 
oxygen alone, we can use the procedure described in Section 3.1 to obtain the corresponding 
values of g(Wmin)- 

Method 1 

This method makes use of the fact that v^^^ > v^l for the same recoil energy, because cal- 

mill rnm '-"■^ ' 

cium is the heavier nucleus and we are only interested in the case where rriy. < mQ. We know 

that g{vmm) must decrease as Vmm increases, so div^i^ i^R)) < di'^^mili^R))- Substituting 
this expression into Eq. (A. 4), it is possible to place a lower bound on the contribution of 
oxygen 

/^y°) CUz(0),A(0))-F(0){Enr /^V^^) ^ .(o,Ca).. ^ /^V^^) 



xtot ■ 



If the differential event rate (dR/dEYi) ° is known, we can make use of this inequality 
and Eq. (A. 3) to infer that 

di?y") fdRY' S^^^'^^HEn) 

dE^J \dEKj l + S(o,c^){En) ■ ^ ' 

In practice, additional factors like the target composition and the acceptance for each nu- 
clear recoil need to be included in 5^'^'*^^^(£'r). In the case of CRESST-II, the minimum 
contribution of oxygen in the relevant energy region is approximately 20%. 

Note, that this method only constrains, but does not actually determine, the ratio of 
oxygen to calcium scatters in each bin. The lower bound that we obtain is almost independent 
of energy although we expect the oxygen fraction to approach 100% as the recoil energy 
increases. To reflect our ignorance of the true ratio, we do not show central values for the 
measurements of g{vrain) obtained with this method. The second method, which we present 
below, will allow for a more accurate determination of the ratio of oxygen to calcium scatters. 

Method 2 

In the approach discussed above, we considered each bin separately and treated the presence 
of a second (heavier) element in the target as an additional uncertainty which reduces the 
information that we can extract from the data on the velocity integral. It would be desirable 
to develop a method that combines the information from different bins in order to use both 
target elements to and constrain the velocity integral more tightly. We now present such a 
method. 

First of all, we must construct energy bins of different sizes according to the following 
rule: Starting from an arbitrary bin [ii^™'"^, i?™*^^] , we use the heavier target (assumed to 
be calcium) to convert the boundaries into t'min-space and then the lighter target (oxygen) 
to convert the velocities back into energy space. This way, we obtain a new energy bin 
E^^\v'^^\Ef'')),E^^\v'^-'^{Ef"'''))\. The two bins thus constructed have the property 
that calcium recoils with energies in bin 1 probe the same region of timm-space as oxygen 
recoils with energies in bin 2. We can repeat this procedure recursively to obtain additional 
bins with increasing bin size (see Figure 14). It is always possible to choose the first bin in 
such a way that the different bins do not overlap. 
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Figure 14. An illustration of the method used to construct bins for CRESST-II. The two bins 
constructed in this way have the property that calcium recoils with energies in bin 1 probe the 
same region of Wnii„-space as oxygen recoils with energies in bin 2. For this plot, we have assumed 
= 9 GeV. 



We only wish to consider bins where a significant excess of signal over expected back- 
ground is observed, meaning that we abort the iteration as soon as we obtain a bin where the 
signal is compatible with the background estimate. To give a concrete example, for niy. = 9 
GeV the energy conversion factor between oxygen and calcium is 

^^ = 1.50, (A.7) 

SO a possible construction for an experiment like CRESST-II would be [8 keV, 12 keV], 
[12 keV, 18 keV], [18 keV, 27 keV]. Unfortunately, the current CRESST-II data reaches 
down only to 10 keV and we cannot bin the data arbitrarily so we will use only the two 
bins [10 keV, 14 keV] and [15 keV, 21 keV] here. For different values of rriy. the conversion 
factor will also change, so different bins must be considered. We present possible choices for 
different DM masses in Table 1. 

We start by evaluating the highest bin and then work our way back towards bin 1. 
In the highest bin, we can reasonably assume that the observed event rate is completely 
dominated by the lightest element in the target. The reason is that we expect events in this 
bin to arise from DM particles close to the escape velocity (otherwise, there would also be 
events in higher energy bins). Consequently, in the corresponding region of fmin-space, we 
expect the velocity integral to decrease rapidly with increasing f mm (see Figure 4) so that a 
lighter target element should have a much larger event rate than a heavier one (probing the 
velocity integral at higher velocities). 

This assumption is good under two conditions. First, that the different target nuclei 
are not too close in mass so that they actually probe different regions of Umin-space and 
second, that there is no additional source of background at large energies that can hide a 
potential signal in higher energy bins. Both of these conditions are fulfilled for the CRESST-II 
experiment, where the background at large energies is approximately flat. 

Under the assumption that only the lightest element contributes in the highest bin, we 
can calculate g(wmin) for this bin in the same way as described in Section 3.1 for a target 
consisting of a single element. Now we can exploit the fact that we have constructed the bins 
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I GeV Conversion factor Bin 1 / keV Bin 2 / keV 

6 L72 10 - 14 17 - 24 

9 1.50 10 - 14 15 - 21 

12 1.34 12 - 15 16 - 20 
15 1.22 



Table 1. The energy conversion factors between oxygen recoils and calcium recoils for constant Umin 
for different choices of the DM mass m^. The third and fourth column give possible choices for 
two bins that are related by this conversion factor. For = 15 GeV, a different binning of the 
CRESST-II data would be required to construct non-overlapping bins. 



in a special way and use the inferred value of ^(fmin) to "predict the number of calcium events 
in the next lower bin. 

We can then subtract the predicted number of calcium events in the next lower bin to 
obtain the number of events that is due to oxygen scatters alone. Again, we are left with 
a number of events that correspond to scattering of the lightest element alone. We can use 
this number to calculate ^(fmin) in the region of Vmin-space corresponding to the second to 
highest bin. Using this measurement to predict the calcium scatters in the next lower bin, we 
can recursively extract g(wmin) for all bins. Compared to the previous method, the advantage 
is clearly a much more accurate prediction of ^(umm)- The disadvantage is that we are forced 
to bin the data in a certain way, which is inconvenient if the energy of each event is not 
publicly available. 



B Constraining the modulation fraction using the absolute rate 

We will demonstrate that it is possible to constrain the modulation fraction once we have 
determined the average value of the rescaled velocity integral from a direct detection exper- 
iment, without assuming any specific properties of the DM halo. The basic idea is that we 
can place an upper bound on the maximum value, ^(vminj summer), and a lower bound on 
the minimum value, ^(fmirn winter), using only the time average, ^(fmin)- 

How such bounds can be obtained is sketched in Figure 15. We observe that 

^(umm, summer) < giy^m - u) (B.l) 
winter) > g{vmm + u) (B.2) 

where u = 29.8 km/s is the velocity of the Earth around the Sun. From this it follows that 

Ag{vtnin) = givmin, Summer) - ^(t^min, winter) < ^(fmin -u) - g{vrmn + u), (B.3) 

and therefore, that the modulation fraction is bounded by 

^ 5 (t^min, summer) - ^(t-min, winter) ^ g{Vmm - u) - g{Vrain + u) 

5(t)mm, summer) +fl'(i'min, wmter) 2g(yrain) 

where we have used giv^m-, summer) + g{vyam-, winter) 2g{ymia)- 

Note that an even stronger constraint can be obtained under the assumption of an 
isotropic DM halo. In this case, only the seasonal change in the radial component of u, Ur, 
matters, because the halo will look the same for any value of and Ud. In this case, we can 
replace u by n,. = ■UC0S7 = 15.0 km/s in the equation above, where 7 is the inclination angle 
between the galactic and ecliptic plane. 
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Figure 15. A graphical representation of the velocity integral. The red shaded region is the region of 
velocity space that is being integrated over, while the white regions are excluded by the requirement 
V > Vmin- We can place an upper bound on the maximum value, g(wmin, summer), and a lower bound 
on the minimum value, g(winin, winter), using only the time average, g(wmin)- Note that the relative 
magnitude of the different velocities are not drawn to scale. 



C Various descriptions of the DM halo 

We present the velocity distributions considered in Sections 4 and 5 and provide the parameter 
choices used for Figures 4 and 8 — for further details see the references. First, we discuss 
isotropic velocity distributions, where f{v) depends only on \v\ and there is no preferred 
direction for DM velocities at the position of the Earth. Then we consider models where the 
anisotropy parameter 

(vo) + {vl) 

/3M.i-^4(^ (CD 

is non-zero, motivated by numerical simulations which favour radially biased orbits with 
/3 = 0.1 — 0.4 at the position of the Earth. 

C.l Isotropic models 
Standard Halo Model 

In the Standard Halo Model (SHM), the velocity distribution function is given by 

f(v) = I ^ [exp(-t'V^o) - ^M-Vesc/^o)] H < Vesc 2) 

where is a normalisation constant. Possible values for the velocity dispersion vq range 
from 180 km/s to 280 km/s, while Vesc can vary between 450 km/s and 650 km/s (see also 
the left panel of Figure 4). Unless explicitly stated otherwise, we will always assume the 
standard values vq = 220 km/s and Vesc = 544 km/s. 

King models 

The cut-off at fesc is introduced by hand in the SHM, which would otherwise predict par- 
ticles with infinitely high velocities and an infinite mass for the Galaxy. This problem is 
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addressed in the King model, which provides a cut-ofF in a self-consistent manner. The 
velocity distribution is obtained from the distribution function: 



where E = '^{x) — jl and ^'(a;) is the relative gravitational potential. The local escape 
velocity at a position x is given by Uesc = A/2^'(a;). At the position of the Earth can 
be parameterised in the same way as the SHM but the parameter a is not directly linked to 
the velocity dispersion and can therefore take values that are much larger than in the SHM 
[46, 47]. Nevertheless, due to its similarity to the SHM we will not discuss King models 
further. 

Double power-lav^ profiles 

A simple modification of the SHM was introduced in Ref.[48]. For double power-law den- 
sity profiles such as the NFW- profile [49], the following ansatz for the velocity distribution 
reproduces better the behaviour at high velocities: 



1 



f{v) = I N [^XP y^tvl ) ' l *^! ^ ^^^'^c . (C.4) 

I 0, > Vc 



\V\ < Ve 



"CSC 



Setting the power-law index k equal to 1 recovers the SHM. The choice 1.5 < k < 3.5 is 
found to give a better fit to velocity distributions extracted from A^-body simulations. We 
use k = 2.5 throughout. 

Tsallis model 

It was argued [34] that the velocity distribution of dark matter particles in numerical simu- 
lations including baryons can be well described by 



1 



l-(l-a)'^ 



1/(1-?) 



f{v) = ^ N ' I'^l < "^esc ^ (c_5) 

see also Ref.[50]. We adopt the parameters q = 0.773, vq = 267.2 km/s and Vesc = 560.8 
km/s from Ref. [34]. 

C.2 Anisotropic models 
Numerical simulations 

A simple anisotropic model has been proposed [5] to describe the data from numerical A-body 
simulations such as Via Lactea [51, 52], GHALO [53] or Aquarius [54]: 

f{v) = I ^ [""""P {-{v^/^rT'-) exp {-{Vt/Vtr')] , 1^1 < ^esc ^ ^^ ^^ 

1 0, [''^1 ^ ''^esc 



where vt = + v"^. For the figures we take the best-fit parameters for the Via Lactea II 

simulation, namely Vr = 202.4 km/s, vt = 128.9 km/s, = 0.934 and at = 0.642 [5], but we 
also show the velocity integral and modulation fraction observed in the GHALOg simulation. 
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Logarithmic ellipsoidal model 

The simplest triaxial generalisation of the velocity distributions considered above was dis- 
cussed in Refs.[2, 55, 56]. We allow a different velocity dispersion in all three directions, 
giving 



1 



f{v) = { ^ f"""^ {-^Ihl - vl/vl - vl/vi^^^ , |^;| < Vesc _ 
i 0, \v\ > Vesc 

The three parameters Vr, v^f, and Vz depend on two constants p and q that describe the density 
distribution and the isotropy parameter 7 (as well as v^). Following Ref.[2], we take p = 0.9, 
q = 0.8 and 7 = —1.33 and calculate Vr, and Vz under the assumption that the Earth is 
on the major axis. 

Distribution functions with /3 = 0.5 

For a constant anisotropy of /3 = 0.5, it is possible to calculate the velocity distribution from 
a given density profile. This was done [57], for centrally cusped density profiles of the form 

For example, for 6 = 4, corresponding to the Hernquist density profile [58], the velocity profile 
is given by 



1 1 _ vl+vl+vl \'^ I I 

^.e.cV^V -c^sc ) ' 1^1 <^-- . (C.9) 

0, \v\ > Vesc 



For 6 = 3, corresponding to the NFW profile [49], an analytical expression of the velocity 
distribution does not exist, but it is straightforward to numerically calculate the velocity 
distribution from the density profile. For the plots shown we adopt a = 10 kpc. 



D Overview of direct detection experiments 

We will briefly discuss the direct detection experiments that we have considered in this 
paper and state the assumptions that we have made. To derive constraints from XENONIOO, 
XENONIO, CDMS-Ge and CRESST-II, we employ the 'maximum gap' method [59], while for 
the low threshold analysis of CDMS-Ge and CDMS-Si, we use the 'binned Poisson' method 
[60, 61]. Our best-fit regions are calculated using a parameter estimation method.^ 

XENONIOO 

We use the most recent data from 100.9 live days of data taking [17]. For the relative 
scintillation efficiency C^s we use a logarithmical extrapolation to zero and calculate the 
energy resolution under the assumption that it is dominated by Poisson fluctuations in the 
number of photoelectrons. 

'For CRESST-II we use Eq. (13) from Ref.[3] to calculate the best-fit parameter region. 
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Bin / keVee 


Event rate 


Error 






in cpd kg^^ keVee~^ 


0.50 - 0.61 


12.8 


1.0 


61 - 72 


10.1 


0.9 


72 - 85 


6.7 


0.8 


U.oO i.Uu 


0.0 


U.o 


1.00 - 1.15 


4.5 


0.7 


Bin / keVee 


Modulation amplitude 


Error 




in cpd kg~^ keVee"^ 


in cpd kg~^ keVee"^ 


0.5 - 1.0 


0.75 


0.54 


1.0 - 2.0 


0.51 


0.30 


2.0 - 3.2 


0.26 


0.16 



Table 2. The CoGeNT event rate and its modulation: the modulation amplitude has been obtained 
by fitting a cosine with period fixed to 365 days and the phase fixed to to = 146 days. 

XENONIO 

It is important to note that the published XENONIO data from Ref.[62] was based on out- 
dated assumptions concerning the relative scintillation efficiency, so all energy scales must 
be rescaled to the same Ces as for the XENONIOO data. We calculate the energy resolution 
in the same way as for XENONIOO. 

We also consider the S'2-only analysis presented in Ref . [63] . We adopt their conservative 
assumption for the ionisation yield with a sharp cut-off at 1.4 keV and take a flat detector 
acceptance of 0.41. To calculate the energy resolution we assume that the production of elec- 
trons is governed by Poisson statistics. We also apply the z-cut suggested by the XENONIO 
collaboration to reduce the number of observed events. 

CoGeNT 

We analyse the publicly available data from the CoGeNT experiment over the course of 
1.2 years [15]. To determine the total event rate, we subtract the L-shell EC contribution 
and a constant background. We use the Lindhard quenching factor and the (corrected) 
detector resolution from Ref. [69]. The efficiency curve was provided by J. Collar (private 
communication). Our bin width is 0.05 keVee. 

To determine the modulation, we fix the peak date to to = 146 days, which is the best- 
fit value for the DAMA modulation. The bins considered and the resulting event rates and 
modulation amplitudes are given in Table 2. 

CDMS-Ge 

We use the final data [18] corresponding to an exposure of 612 kg days. We assume a detector 

resolution of 0.4 keV and take the detector acceptance from Ref. [64]. 

We also consider the dedicated low-threshold analysis of the CDMS-Ge data [19] and use 

only the data of the most constraining detector, T1Z5. Here, we adopt an energy resolution 
of 0.2(SR-keV)V2 [g5]^ 

CDMS-Si 

To obtain a constraint from the silicon detectors in CDMS-II, we take the low-threshold data 
from the shallow-site run [66] . We include the runs at 3V and at 6V but consider only the Z4 
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Bin / keV 


Total events 


Expected background 


Signal (min) 


Signal (max) 


10 - 13 


9 


3.2 


3.2 


9.6 


13 - 16 


15 


6.1 


5.3 


13.3 


16 - 19 


11 


7.0 


1.2 


7.8 


19 - 22 


8 


6.3 


0.2 


5.0 


22 ~ 25 


4 


5.2 





1.8 


25 - 28 


4 


4.6 





2.2 


28 - 31 


4 


4.3 





2.5 


31 - 34 


3 


4.0 





1.5 


34 - 37 


4 


3.7 





3.0 


37 - 40 


5 


3.5 


0.2 


4.3 


13 - 15 


11 


3.8 


4.0 


11.0 


16 - 19 


11 


7.0 


1.2 


7.8 


19 - 22 


8 


6.3 


0.2 


5.0 


10 - 14 


18 


5.4 


8.4 


17.4 


12 - 15 


19 


5.7 


9.1 


18.1 


15 ~ 21 


19 


13.7 


1.7 


10.1 


16 - 20 


13 


9.4 


0.9 


7.9 


17 - 24 


17 


14.1 


0.5 


7.7 



Table 3. Data from the CRESST-II experiment: total mimber of events, expected number of back- 
ground events and la confidence interval for the signal using a Feldman-Cousins approach. 

detector which gives the strongest bound. Resolution and detector efficiency are both taken 
from Ref. [66]. 

SIMPLE 

We do not attempt to combine the results from Stage 1 [67] and Stage 2 [68] but consider 
only the latter. Since the recoil energy is not measured, the experiment can place only an 
upper limit on the total event rate above the detector threshold of 8 keV. We take the cut 
acceptance and nucleation efficiency from Ref. [68]. 

CRESST-II 

We use the recent results from 730 kg days of data taking [16]. We estimate that the 
acceptance region corresponds to 86% acceptance for oxygen, 89% for tungsten and 90% for 
calcium. We take the energy resolution to be 0.3 keV for each detector. To extract the DM 
signal, we read off the total number of measured events and the expected background from 
Figure 11 of Ref. [16] and use the Feldman-Cousins approach [70]. 

Depending on the application, we bin the CRESST-II data in different ways. In order to 
calculate best-fit parameter regions, we consider a bin width of 3 keV, rather than a bin width 
of 1 keV as presented in Ref. [16]. However, if we want to calculate ^(t'min) using Method 1 
from Appendix A. 2, we have to make sure that the detector efficiency does not vary strongly 
within a single bin due to the onset of additional detector modules. This condition is fulfilled 
for the three bins [13 keV - 15 keV], [16 keV - 19 keV] and [19 keV - 22 keV]. Finally, to 
apply Method 2 from Appendix A. 2, we need to bin the data in a special way, namely in the 
two bins [10 keV - 14 keV] and [15 keV - 21 keV]. The resulting number of events is shown 
in Table 3. 
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Finally, we also consider the improved constraint from the CRESST-II commissioning 
run presented in Ref. [71] using the detector efficiencies given there. 

DAMA 

We use the combined observed event rates from DAMA/Nal and DAMA/LIBRA [14] for the 
8 bins spanning the energy range 2-6 keVee. No significant modulation is observed at higher 
energies. As we are interested in DM particles with mass below 20 GeV, we neglect scattering 
on iodine and consider only the contribution of sodium. We assume a quenching factor of 
0.3 (other values are discussed in Refs.[24, 45]) and no contribution from channelling [72]. 
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